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ABSTRACT 


Based  on  the  three-dimensional  theory  developed  at  Therm  Advanced 
Research,  a  simple  procedure  for  calculating  the  detailed  aerodynamic 
loading  on  the  shroud  of  a  finite-bladed  ducted  propeller  in  forward 
flight  is  presented. 

This  procedure  has  been  designed  for  the  evaluation  of  the  several 
thousand  configurations  exactly  represented  by  the  tabulated  values  of 
the  pertinent  parameters,  or  reasonably  approximated  by  them.  In  ad¬ 
dition,  configurations  whose  parameters  fall  somewhat  in  between  or 
outside  these  values  can  be  evaluated  by  suitable  interpolation  or  ex¬ 
trapolation.  The  data  were  selected  through  liaison  with  propeller  and 
aircraft  manufacturers  to  encompass  the  current  state  of  the  art  for 
ducted  propeller  design. 

As  a  result,  it  is  possible  to  carry  out  a  number  of  calculations  by 
hand  quite  readily  which  could  not  be  done  heretofore,  for  example,  the 
examination  of  the  effect  of  the  propeller  advance  ratio. 
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parameter  for  NACA  4-digit  mean  line,  y  /c  ; 
equivalent  to  m  of  Ref.  13  ^ 


shroud  sectional  pitching  moment  about  quarter 
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shroud  thickness 


Legendre  function  of  the  second  kind  and 

half  integer  order  with  argument  co  ,  see  Ref.  14 


shroud  sectional  radial  force  or  radial  force  per 
unit  circumferential  length,  positive  outward 


shroud  reference  radius,  taken  as  radius  of 
camber  line  at  propeller  plane,  see  WORKSHEET  I 


radius  of  shroud  inner  surface  at  propeller  plane, 
see  WORKSHEET  I 


propeller  radius,  see  WORKSHEET  I 


elements  for  shroud  curvature  contribution  to 
shroud  surface  pressure  arising  from  effective 
shroud  camber 


total  axial  force  on  shroud,  positive  upstream 


maximum  shroud  thickness,  see  WORKSHEET  I 


shroud  thickness  at  propeller  plane,  see 
WORKSHEET  I 


propeller  thrust 


elements  for  sliroud  curvature  contribution  to  shroud 
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axial  position  of  aft  end  of  constant  portion  of 
two-dimensional  loading  obtained  from  NACA  6-digit 
mean  lines;  measured  relative  to  midchord,  see 
WORKSHEET  I 


axial  position  of  shroud  sectional  center  of  pres¬ 
sure;  measured  relative  to  midchord,  see 
WORKSHEET  I 


axial  position  of  maximum  ordinate  of  NACA  4-digit 
mean  lines  measured  along  shroud  chord  line;  here 
taken  relative  to  midchord  in  linearized  sense,  see 
WORKSHEET  I 


axial  position  of  propeller  plane;  measured  relative 
to  midchord,  see  WORKSHEET  I 
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viii 
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propeller  circulation  distribution 
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e(x/c) 


two-dimensional  net  pressure  coefficient  for  NACA 
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chord,  X  /c 
’  P 


angular  velocity  of  propeller 


1 

;j 


A  DETAILED  NUMERICAL  EVALUATION 
OF  SHROUD  PERFORMANCE 
FOR  FINITE-BLADED  DUCTED  PROPELLERS 

Introduction 

Under  the  sponsorship  of  the  Office  of  Naval  Research,  personnel 
of  Therm  Advanced  Research  have  been  engaged  for  the  past  four  and 
one-half  years  in  the  general  study  of  ducted  propellers  with  finite 
blade  number.  As  part  of  a  long  range  ONR  V/STOL  program,  the 
basic  objective  has  been  to  provide  a  definitive  analysis  for  different 
flight  regimes. 

A  theory  of  the  forward  flight  regime  for  zero  angle  of  attack  has 
1  “6 

been  completed  ,  including  the  effects  of  shroud  camber,  thickness 
and  chord  to  diameter  ratio,  as  well  as  the  propeller  loading,  blade 
number,  advance  ratio,  tip  clearance  and  axial  position.  These 
results  furnish  not  only  general  formulas  for  both  the  steady  and 
harmonic  shroud  loads  but  also  considerable  physical  insight  into 
the  overall  problem.  However,  the  analytic  solution  is  not  amenable 
to  detailed  aerodynamic  evaluation  because  it  requires  the  use  of 
electronic  computers  to  perform  certain  numerical  calculations. 

Since  the  engineering  value  of  any  such  solution  is  measured  by  the 
ease  with  which  data  can  be  extracted,  these  machine  calculations  were 
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carried  out  at  TAR  for  representative  ducted  propeller  configurations. 
This  report  is  the  result:  namely,  the  development  of  a  simple  pro¬ 
cedure  whereby  the  steady  shroud  net  loading,  surface  pressures, 
thrust,  sectional  radial  force,  sectional  pitching  moment  and  sectional 
center  of  pressure  can  be  calculated  by  hand. 

To  review  briefly,  the  theory  was  derived  as  follows.  The  fluid  is 
taken  as  inviscid  and  incompressible.  Viewed  in  propeller-fixed  coordi¬ 
nates,  the  duct  is  replaced  by  suitable  distributions  of  sources  and 
vortices  on  a  cylindrical  reference  surface  whose  axis  is  coincident  with 
the  duct  axis.  The  propeller  blades  are  represented  by  individual, 
radial  vortex  lines  of  varying  circulation  and  accompanying  trailing 
helical  vortex  sheets.  K  the  form  of  the  propeller  circulation  is  as¬ 
sumed,  the  formulation  of  the  governing  equations  then  follows  the  pro¬ 
cedure  of  thin  lifting- surface  theory.  Within  the  limitations  of  this 
theory,  we  find  that  the  problem  for  the  steady  load  is  equivalent  to  an 

7 

axisymmetric  ring  wing  with  the  same  thickness  distribution  as  the 
given  duct  but  an  "effective”  camber  composed  of  the  original  geometric 
camber  plus  the  radial  wash  of  the  shroud  sources  and  propeller  trailing 
vortices. 

The  source  strength  is  simply  proportional  to  the  rate  of  change  of 

airfoil  thickness  and  its  effects  are  quite  readily  adaptable  to  a  hand 

5  6 

calculation  procedure  ’  .  To  determine  the  strength  of  the  shroud 


vortices  on  the  other  hand,  it  is  necessary  to  calculate  the  two- 
dimensional  Glauert  coefficients  of  the  effective  camber  and  correct 
them  in  turn  for  shroud  curvature.  The  curvature  correction  is  no 
difficulty  but  the  calculation  of  the  Glauert  coefficients  is,  i.  e. ,  a 
technique  for  numerical  Fourier  analysis  is  required.  Unfortunately 
the  nature  of  the  propeller  contribution  to  the  effective  camber  in  par¬ 
ticular  is  such  that  the  available  hand  methods,  e.g. ,  Ref.  8,  are  in¬ 
adequate.  As  a  result,  the  Fourier  analysis  had  to  be  carried  out  by 
machine  and  this,  together  with  the  propeller  camber  itself,  constituted 
the  major  computational  effort. 

With  regard  to  the  propeller  circulation,  distributions  for  a  so- 
called  optimum  without  a  hub  are  used.  NACA  thickness  and  camber 
distributions  have  been  chosen  for  the  shroud  airfoil  section.  Data  are 
tabulated  for  several  values  of  the  pertinent  parameters  which  encom¬ 
pass  the  current  state  of  the  art  for  ducted  propeller  design.  Utilizing 
all  possible  combinations  of  these  tables  enables  the  exact  evaluation  of 
several  thousand  ducted  propeller  configurations.  Reasonable  "non¬ 
exact"  caser  '•an  be  evaluated  by  proper  interpolation  or  extrapolation. 

Chapter  One  elaborates  on  the  choice  and  other  details  of  the 
propeller  circulation,  the  shroud  airfoil  sections,  and  the  range  of  the 
gross  parameters  of  geometry  and  flight  condition.  The  method  of 
computation  and  related  information  for  the  tables,  the  format  of  the 


evaluation  and  a  suggested  interpolation  or  extrapolation  technique  are 
also  presented.  Chapter  Two  contains  illustrative  numerical  examples 
giving  all  the  work  as  it  should  appear  on  the  worksheets.  The  report 
concludes  with  the  actual  evaluation  procedure  and  the  tabulated  data. 

It  is  hoped  that  this  design  procedure  will  pave  the  way  for  de¬ 
tailed  evaluation  and  parametric  study  of  ducted  propeller  configura¬ 
tions  which  were  impracticable  previously.  It  is  also  hoped  that  this 
report  will  be  recognized  as  a  source  of  information  of  a  more  universal 
nature,  e.  g. ,  the  Glauert  coefficients  for  the  NACA  4-,  5-  and  6-digit 
mean  lines. 


CHAPTER  ONE 


GENERAL  DISCUSSION 

Propeller  Circulation  Distributions 

The  form  of  the  distribution  of  the  propeller  circulation  r  deter¬ 
mines  the  propeller  induced  camber  and  the  direct  contribution  of 
the  propeller  c  to  the  shroud  surface  pressure  coefficient  c 

r  p 

There  are  several  choices.  As  opposed  to  our  original  thought^,  we 

g 

have  picked  the  classical  Betz  optimum  with  a  tip  correction  derived 
by  T.  Goodman^*^  to  account  for  the  proximity  of  the  duct.  Since  cer¬ 
tain  three-dimensional  effects  of  both  the  propeller  and  shroud  are 
omitted,  such  an  optimum  will  not  be  a  true  one.  Further  difficulties 
may  arise  from  viscous  effects  if  the  tip  lies  within  the  shroud  boundary 
layer.  Despite  these  complications,  however,  it  should  yield  fairly 
accurate  results  for  the  propeller- shroud  interaction  phenomena  and 
moreover,  it  is  relatively  easy  to  compute,  cf.  Ref.  11. 

The  specific  circulation  distributions  used  in  the  evaluation  proce¬ 
dure  are  given  in  TABLES  1.1  -  1.3,  each  normalized  so  as  to  produce  a 
unit  propeller  thrust  coefficient  C.p  =  1  neglecting  inflow.  The  other 
parametric  inputs  consist  of  the  number  of  blades  N  ,  the  ratio  of  the 
propeller  radius  to  the  radius  of  the  inner  shroud  surface  e  =  > 

and  the  propeller  advance  ratio  or  ratio  of  the  forward  speed  to  the  tip 
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speed  J  =  U/S2Rp  .  Because  of  the  square-root  behavior  of  these 
distributions  at  the  tip,  a  Glauert-type  angular  variable  /S  was  intro¬ 
duced  to  insure  sufficient  accuracy.  Though  increments  of  only  5°  are 
shown,  increments  down  to  0. 125°  were  computed  in  some  cases.  The 
corresponding  radial  stations  r/R  are  given  for  convenience. 

These  circulation  distributions  are  not  used  directly  in  the  proce¬ 
dure  but  are  presented  for  general  information  and  completeness.  For 
example,  if  evaluation  for  a  different  bladewise  thrust  gradient 

r(r/Rp) 

dr/Rp  RpU 

is  desired,  it  can  be  approximated  by  superposition  though  this  capa¬ 
bility  has  not  been  incorporated  explicitly  in  the  procedure. 

When  a  hub  is  present,  the  shape  of  the  distribution  is  altered 
such  that  it  cannot  be  fitted  by  the  given  distributions.  However,  these 

changes  are  relatively  small  and  essentially  confined  to  a  region  near 
12 

the  hub  .  To  simulate  a  hub  in  a  crude  fashion,  configurations  have 
been  run  in  which  the  tabulated  circulation  distributions  have  been 
simply  truncated  at  the  hub.  The  results  indicate  that  hubs  with 
diameters  up  to  about  257oOf  the  propeller  diameter  have  a  negligible 
effect  on  the  shroud  performance  and  that  the  performance  can  be 


evaluated  using  the  given  circulation  distributions.  The  direct  effect 
of  the  radial  velocity  induced  by  the  hub  depends  upon  the  shape  and  ex¬ 
tent  of  the  hub  and  can  be  estimated  from  slender  body  theory. 

The  results  for  the  case  of  a  uniform  circulation  have  been  com¬ 
pared  with  the  tabulated  distributions.  Since  they  are  independent  of 
e  and  for  a  fixed  value  of  C.p  ,  also  N  and  J  ,  the  number  of  calcu¬ 
lations  are  greatly  reduced  but  some  discrepancies  are  introduced.  In 
general,  the  two-dimensional  Glauert  coefficients  of  the  propeller  in¬ 
duced  camber  e_ ,  fall  within  the  tabulated  values  for  a  given  set  of 
overall  geometric  parameters  to  be  described  later.  The  discrepancy 
in  the  values  averages  about  5% .  For  extremum  values  of  J  ,  though, 
the  discrepancy  can  be  as  large  as  20% .  The  same  magnitude  of  dis¬ 
crepancy  is  also  introduced  in  the  direct  propeller  contribution  to  the 
shroud  surface  pressure  coefficient  c^^  ,  and  the  propeller- shroud 
thickness  contribution  to  the  shroud  thrust  coefficient  ^  . 

Shroud  Thickness  and  Camber  Distributions 

In  the  evaluation  procedure  typical  shroud  thickness  and  camber 

distributions  are  represented  by  the  NACA  4-,  5-  and  6-digit  series  of 

airfoil  sections  and  mean  lines.  These  distributions  have  been  selected 
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because  they  are  well  known  and  data  for  them  is  readily  available  . 
Other  shapes  may  be  obtained  by  superposition.  Data  is  also  included 
for  the  idealized  case  of  zero  thickness. 


The  scaling  rules  appropriate  to  each  NACA  airfoil  and  mean  line 
have  been  incorporated  in  the  design  procedure.  However,  for  6-digit 
airfoil  sections  which  can  be  scaled  only  within  a  given  family  and  for 
5-digit  mean  lines  which  can  not  be  scaled,  the  evaluation  procedure 
contains  only  data  for  those  distributions  which  are  listed  in  the  TABLE 
OF  PARAMETERS. 

In  general,  standard  NACA  symbols  are  preserved,  but  those 
parameters  which  depend  upon  chordwise  position  are  referred  to  the 
midchord  of  the  shroud  instead  of  the  leading  edge  in  order  to  be  consis¬ 
tent  with  the  coordinate  system  we  have  adopted,  see  WORKSHEET  I. 
The  4-  and  5-digit  airfoil  sections  are  characterized  by  the  ratio  of  the 
maximum  thickness  to  the  chord  t^^/c.  The  6-digit  airfoil  sections  are 
characterized  by  their  family,  i.  e. ,  the  first  two  digits  of  the  NACA 
designation;  within  a  given  family,  the  airfoil  is  characterized  by  the 
ratio  t^/c.  The  4-digit  mean  lines  are  characterized  by  the  ratio  of 

the  maximum  mean  line  ordinate  to  the  chord  m  =  y  /c  and  the  ratio 

•'m 

of  the  axial  position  of  this  ordinate  to  the  chord  p  =  x^^/c.  The  5- 
digit  mean  line  is  characterized  by  its  NACA  designation.  The  6-digit 
mean  line  is  characterized  by  the  ratio  of  the  axial  position  of  the  aft 
end  of  the  constant  portion  of  the  two-dimensional  loading  to  the  chord 

=  x^/c  and  the  two-dimensional  design  lift  coefficient  c^  . 

i 
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3  Gross  Geometry  and  Flight  Parameters 

As  we  have  seen,  the  basic  component  parameters  for  the  propeller 
circulation  are  C,p  ,  N  ,  e  and  J  ;  for  the  shroud  thickness  distribu- 
tion,  t^/c  or  the  NAG  A  family  designation;  and  for  the  shroud  geo¬ 
metric  camber  distribution,  m  and  p  ,  or  the  NACA  designation,  or 

a,  and  c^  .In  addition,  there  are  three  other  parameters  for  the 
i 

ducted  propeller  as  a  whole;  the  relative  location  of  the  propeller  plane 
with  respect  to  the  midchord  x  =  x^/c  ,  the  shroud  chord  to  diameter 
ratio  X  =  c/2R  ,  where  R  is  the  reference  radius  of  the  shroud 
camber  line  at  the  propeller  plane,  and  the  ratio  of  the  propeller  radius 
to  the  reference  radius  jii  =  thickness  distribution  param- 

sr 

eter  t^/c  is  not  an  independent  parameter  of  the  problem.  Once  the 
values  of  X  ,  e  ,  and  p  are  specified,  the  value  of  t^^/c  for  a 
given  NACA  airfoil  section  is  uniquely  determined.  In  order  to  elimi¬ 
nate  the  possibility  of  using  the  incorrect  value  of  t^/c  ,  the  evalua¬ 
tion  procedure  contains  a  method  whereby  the  proper  value  is  always 
computed. 

The  range  of  the  parameters  chosen  is  shown  in  the  TABLE  OF 
PARAMETERS.  Certain  of  the  parameters  take  on  only  specific  values 
while  others  are  arbitrary.  Data  has  been  generated  and  is  presented 
for  all  possible  combinations  of  the  indicated  values.  In  view  of  the 
large  number  of  parameters,  considerable  selectivity  had  to  be 


exercised.  To  make  the  values  chosen  as  representative  of  ducted  pro¬ 
peller  practice  as  possible,  two  things  were  done.  First  we  made  an 
extensive  compilation  of  the  parameters  for  all  the  known  configurations 
which  have  been  tested  and  on  this  basis  made  our  preliminary  choices. 
This  was  followed  up  by  discussion  and  contact  with  various  propeller 
and  aircraft  manufacturers.  Their  thoughts  then,  together  with  consid¬ 
eration  of  what  was  a  reasonable  computational  effort,  guided  us  to  the 
final  values  adopted. 

Computation  of  Tables 

Most  of  the  tables  were  generated  specifically  for  the  present 
report.  No  derivations  are  given  but  the  equations  used  are  shown  on 
the  cover  sheet  which  precedes  each  set.  If  necessary,  these  equations 
may  be  obtained  from  the  cited  references  generally  in  the  exact  form 
as  shown.  In  the  cases  where  the  form  is  different,  an  obvious  change 
in  variable  and/or  integration  by  parts  is  required.  The  rest  of  the 
tables  were  taken  from  previous  reports  and  are  simply  reproduced. 

Also  included  on  the  coversheet  is  the  accuracy  of  the  numbers  and 
any  pertinent  information  regarding  the  use  of  the  tables.  For  example, 
on  the  cover  sheet  for  TABLES  6. 1  -  6. 18,  there  are  instructions  for 
the  use  of  the  tables  when  x  =  0. 25  .  These  instructions  utilize  certain 
symmetries  in  order  to  economize  on  the  number  of  tables  which  have  to 
be  presented.  It  is  thus  important  that  the  coversheets  be  read  carefully 


before  using  the  tables. 

All  new  machine  computations  were  carried  out  on  the  IBM  7090  or 

CDC  1604.  The  incomplete  elliptic  integral  of  the  first  kind  F(v’  ,k)  was 

computed  by  a  modified  SHARE  program  and  the  Legendre  function  of  the 
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second  kind  and  half  integer  order  ^  o^r  own  program  . 

A  Generalized  Simpson’s  Rule  was  written  for  all  integrals.  In  this  pro¬ 
gram  the  interval  of  integration  is  successively  halved  until  three  con¬ 
secutive  answers  are  obtained  which  agree  to  within  the  prescribed 
accuracy. 

We  have  checked  and  proofed  the  tables  thoroughly.  All  equations 
have  been  double-checked  for  analytical  correctness  by  rederivation,  as 
well  as  all  computer  programs  by  both  spot  and  sample  hand  calculations 
using  numerical  procedures  different  from  the  machine.  Proof-reading 
has  been  done  three  times  to  insure  that  all  data  was  transcribed  ac¬ 
curately  from  the  machine  printout. 

Format  of  Evaluation  Procedure 

The  evaluation  procedure  consists  of  six  (6)  sub-procedures  and 
two  worksheets  on  which  all  calculations  are  done  directly.  The  sub¬ 
procedures  are: 

1.  DETERMINATION  OF  PARAMETERS 

2.  COMPUTATION  OF  THREE-DIMENSIONAL  GLAUERT 
COEFFICIENTS  OF  EFFECTIVE  CAMBER 


3.  COMPUTATION  OF  SHROUD  SECTIONAL  RADIAL 
FORCE  AND  MOMENT  COEFFICIENTS  AND  CENTER 
OF  PRESSURE 

4.  COMPUTATION  OF  SHROUD  THRUST  COEFFICIENT 

5.  COMPUTATION  OF  NET  SHROUD  PRESSURE 
COEFFICIENTS 

6.  COMPUTATION  OF  OUTER  AND  INNER  SHROUD  SURFACE 
PRESSURE  COEFFICIENTS 

Two  sample  prints  of  each  worksheet  accompany  the  report 
for  trial  calculations.  Like  the  tables,  the  previous  references 
should  be  consulted  for  the  basic  background  on  the  manipulations 
involved. 

The  order  in  which  the  computations  are  done  is  essentially 
arbitrary,  except  for  SUB -PROCEDURES  1  and  2.  These  must  always 
be  done  first  and  in  sequential  order  regardless  of  what  other  results 
are  required.  Once  SUB-PROCEDURE  1  is  completed,  any  of  the 
remaining  may  be  done  by  non- engineering  personnel. 

The  sub-procedural  instructions  indicate  what  tables  should  be 
taken  and  where  to  write  them  on  the  worksheets.  The  worksheets  are 
made  up  of  an  appropriate  layout  of  the  ducted  propeller  geometry  and 
numbered  columns.  The  worksheets,  themselves,  also  contain  a  set 
of  instructions  which  indicate  the  mathematical  operations  to  be  per¬ 
formed  after  completion  of  the  respective  sub-procedural  steps.  These 
instructions  are  written  primarily  in  equation  form  using  the  column 


headings  as  the  terms  of  the  equations.  The  indicated  operations  are 
to  be  performed  on  each  element  in  the  column.  The  resulting  directions 
of  action  of  all  forces  and  moments  are  shown  pictorially. 

1. 6  Interpolation/Extrapolation 

The  evaluation  procedure  has  been  designed  for  the  evaluation  of 
the  configurations  exactly  represented  by  the  particular  values  chosen 
for  the  parameters  or  reasonably  approximated  by  these  values. 
Nevertheless,  "non-exact"  cases,  i.  e. ,  configurations  whose  param¬ 
eters  fall  somewhat  in  between  or  outside  those  tabulated,  can  be 
evaluated  by  suitable  interpolation  or  extrapolation. 

Experience  has  shown  us  that  the  best  way  is  first  to  evaluate 
several  "exact"  configurations  using  values  of  the  tabulated  parameters 
which  bound  or  adjoin  those  for  the  "non-exact"  case  and,  then,  to 
interpolate  or  extrapolate  in  the  final  results.  The  actual  value  of  the 
arbitrary  parameters  are  used  in  the  evaluations  with  the  exception  of 
tj^/ c.  The  value  of  tj^/ c  used  must  be  consistent  with  the  values  of 
A,  e  and  [i  as  explained  previously  on  p.  9.  Although  this  technique 
requires  the  evaluation  of  several  configurations,  it  has  one  great 
advantage.  Namely,  all  other  "non-exact"  cases  anywhere  in  the 
same  region  of  interpolation  or  extrapolation  are  readily  found. 

Direct  interpolation  or  extrapolation  of  the  tabulated  data  before 
evaluation  will  also  work,  provided  only  one  of  the  actual  parameters 
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does  not  correspond  to  the  tabulated  parameters.  If  more  than  one 
parameter  is  different,  then  one  or  more  multiple  interpolations  or 
extrapolations  are  generally  required.  This  becomes  not  only  very 
complicated  and  time  consuming,  but  results  in  the  loss  of  accuracy 
through  both  computational  and  arithmetical  errors. 


CHAPTER  TWO 


ILLUSTRATIVE  EXAMPLES 

Background 

Three  typical  examples  are  presented  to  illustrate  some  of  the 
ways  in  which  the  evaluation  procedure  can  be  applied,  the  numerical 
operations  involved  in  each  sub-procedure  and  the  use  of  the  tables. 

For  each  example  a  set  of  completed  worksheets  are  given  showing 
the  proper  number  of  decimal  places  consistent  with  the  accuracy  of 
the  tables.  We  have  kept  three  places  for  the  most  part.  Where 
readily  practicable,  though,  we  have  retained  four  places  to  minimize 
the  accumulation  of  round-off  errors. 

Example  for  "Reasonably  Approximate"  Configurations 

The  first  example  is  the  complete  evaluation  of  a  configuration 
which  is  reasonably  approximated  by  the  tabulated  data.  The  com¬ 
pleted  worksheets  are  reproduced  in  Figs.  1  and  2  and  can  be  worked 
out  straightforwardly  by  following  the  sub-procedures.  The  dimensions 
given  correspond,  of  course,  to  the  actual  configuration,  but  the  dimen¬ 
sions  for  the  "reasonably  approximate"  configuration  are  found  by 
modifications  consistent  with  the  chosen  tabulated  parameters.  For 
the  same  propeller  radius,  the  dimensions  of  this  configuration  are 
Rp  =  3. 5000  ft  ,  c  =  3. 8889  ft  ,  R.  =  3. 5354  ft  , 
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SUB-PROCEDURE  2 
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X  =  -0.9722  ft  5  and  t  =  0.7070  ft  .  The  "new"  value  of  t 
P  P  P 

used  in  the  calculations  was  0.  6614  ft  ;  see  step  4  of  SUB-PROCE¬ 
DURE  1.  It  should  also  be  noted  that  the  value  of  changes  to  140 
rad/sec.  The  respective  configurations  are  drawn  to  scale  in  Fig.  3 
and  do  not  appear  appreciably  different. 

An  ’’exact"  evaluation  of  the  actual  configuration  has  also  been  run 
by  means  of  the  general  computer  programs.  The  worksheets,  com¬ 
pleted  as  if  tables  were  available  for  these  values  of  the  parameters, 
are  shown  in  Figs.  4  and  5.  For  this  case  the  tabulated  data  for 


sheet  for  TABLES  1. 1  -  1. 3  with  F  (<p,k)/F(7r/2,k)  being  obtained  by 
four -point  Lagrangian  interpolation  of  data  not  published  in  this  report 
but  which  is  available  in  the  general  computer  program.  Similarly, 
the  elements  ^  were  obtained  bj'^  three -point  Lagrangian  interpola¬ 
tion  of  data  in  Ref.  2,  and  the  elements  ^  and  Tj^  ^  by  three - 
point  Lagrangian  interpolation  of  the  appropriate  tables  contained  in 
this  report.  Within  the  range  of  the  parameters  N  ,  e  ,  J  and  X  , 
the  interpolated  results  are  as  accurate  as  those  which  can  be  obtained 
by  actual  calculations  using  the  equations  on  the  appropriate  cover- 
sheet. 

The  relative  magnitudes  of  the  shroud  sectional  force  and  moment 
coefficients  and  center  of  pressure,  the  shroud  thrust  coefficient  and 


(See  Figs.  1,  T  4  and  5  for  comparison  of  various  coefficients.) 


7?  79  74  I  “N 

0.99  7  -P.PIO  0.013  j 


Write  the  «t«mcnu  of  couv'iOS. 
corretpondlnx  to  ('‘-2.9.4...  12, 
aa  the  eorrea ponding  elements 
of  cou-MM  86. 

©  Write  the  first  element  of 

rou-MV  93<  eorrespondlng  to 
i'  0  ,  AS  the  element  of 
COWMS  89, 

©  Write  the  seeond  element  of 
couMV  39,  corresponding  to 
S'  I  ,  as  the  element  of 
COU'MM  90. 

Write  the  elements  of  cou'M'iOI, 
corresponding  to  v  2,3.4...  12. 
as  the  corresponding  elements 
Ot  VOU'U'l  91. 


5  as  the  elements  of  coux't  79. 

Write  the  elements  of 
cou-M*  61  •  corresponding 
(oi"^l.2.9...  II,  as  the 
corresponding  elements 
of  eoa-M-e  87. 

Write  the  first  element  of 
cotevs  61,  eorresponding  to 
t"0  a  as  the  element  of 
eoUKS  91. 

Write  the  second  element  of 
couvs  61,  corresponding  to 
t'^l  .  as  the  element  of 
cou-Ms  92. 
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0.0291 
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3 

o.ootl 

4 

0.0091 
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0.0091 
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O.OO^j 

89 

90 

— 2l 

’0.062 

0,039 

o.oJ 

_s 

75  76  77 

•0.741  .0.250  .0.991 


@  80-78x79 

82  *  80  X  81 
85  »  83  X  84 
88  <  86  X  87 


©91 '.89x91 
^  96  •89x92 
98  *  90x92 

100  *  Algdbrale  sum  of  85 


FIGURE  4 

COMPLETED  WORKSHEET  I  FOR  "EXACT"  EVALUATION  O 
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SUB-PROCEDURE  2 

COMPUTATION  OP  THREE-DIMENSIONAL  GLAUERT  COEPnCIENTS  OF  EFFECTIVE  CAMDER 
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3 

4 

5 
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1.4239 
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0.1863 

0.2661 

i 

>3. 1513 

•0.5934 

2 

-1.6626 

-0.3166 

3 

2.6616 

0.5426 

4 

-5.2139 

•0.9616 

5 

0 

0 

6 

0 

0 

© 


Write  the  dement*  of  couvi  &  as  the  first 
element  of  cou-MM  )«  9.  11.  13.  15.  17  and  19 
respectively. 
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20*6X7 

21- 0x0 

22- 10x11 

23  •  12  X  13 

24  -  14  X  IS 

25  •  10  K  17 

20-t0xi» 

34  20  *  21  •  22  «  23  »  24  *  2S  >  26 


I 


10 

0 

•0.3166 

0.0466 

0 

-0.0092 

0 

-0.0013 

0 

12 

13 

0.0257 

0,5426 

0 

•0.0019 

0 

0.0019 

0 
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14 

15 

0 

•0.9616 

0.0053 

0 
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0 
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0 
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0 

0 
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0 
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0 

0.0001 

1# 
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0 

0 

0.0003 

0 

0.0003 

0 
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0 

20 

21 

22 

23 

24 

25 

26 

0.2017 
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0 
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0 

•0.0146 

0 
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0 

0 
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0 
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0 

0 

0 
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0 

0.0029 

0 
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0 
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0 
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0 

0 
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0 

0 

0 
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30 

31 

32 

33 
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35 

36 

37 

36 

39 
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•0.025 
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•0.014 
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0,014 
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7 

-0.003 

0,001 

0.000 

•0.003 

0.000 

0 

•0.003 

«. 

•0.001 

•0,002 

‘0.001 

•0.001 
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0 
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II 

1  0.000 
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0 
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Write  the  first  seven  elements 
of  coU'Mx  39  as  the  first  element  of 
eou^vs  41.  43,  45.  47,  49.  51  and  53 
rcspeelhely. 
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54  •  40  X  41 

55  *  42x43 
56-44X45 
57*  46X47 
56*46X49 
59  •  50  X  SI 
60-52X53 
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Q  cj*®*6l*54*55*56*57*58*59*  60i  t'-0,1.2,„6 
«t.*®-6l*39;  t'-t.O.O...  12 
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4 

O.vO# 

5 
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6 
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10 
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SUB-PROCEDURE  4 

COMPUTATION  OF  SHROUD  THRUST  COEPHCIENT 


pi  as  the  elements  of  coukx  79. 

?  Write  the  elements  of 

'  rou-xt  61.  corresponding 

toP«l,2.3  ...  II,  as  the 
corresponding  elements 
of  couxx  67. 

^  Write  the  first  element  of 
cou-wr  61.  corresponding  Co 
I'^'O  >  as  the  element  of 
cou-wr  91. 

^  Write  the  second  element  of 
COUNTS  61,  corresponding  to 
p  *1  •  as  the  clement  of 
cou-Mt  92. 
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79 

60 

81 
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0.300 
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107  102-  103*  101*  105*  106 
106  101 X  107 

109  '  Algebraic  sum  of  e2 
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|gure  4 

I*”  evaluation  of  actual  configuration 


SUB-PROCEDURE  5 

COMPUTATION  OP  NBT  SHROUD  PRESSURE  COBPnClENTS 


I  Writ*  th«  eltmMU  of  couvN  i\ 

OB  WORKSHEET  1  «•  the  firet  element  of 
cou-MNi  lit.  Itl.  123.  12$,  127.  12t. 
191.  199.  19$.  137.  13».  HI  and  143 
reapcctlvely. 
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j^LUATION  OF  ACTUAL  CONFIGURATION 
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the  shroud  outer  and  inner  surface  pressure  coefficients  for  the  two 
cases  may  be  found  from  Figs.  1  and  4  and  Figs.  2  and  5,  respectively. 
The  overall  numerical  agreement  indicates  roughly  the  level  of  con¬ 
fidence  which  can  be  attached  to  the  use  of  the  evaluation  procedure 
for  "reasonably  approximate"  configurations.  In  general,  it  appears 
chat  we  obtain  a  satisfactory  engineering  estimate. 

2.  3  Example  for  Interpolation/Extrapolation 

The  second  example  is  the  computation  of  the  outer  and  inner 
shroud  surface  distributions  for  a  "non-exact"  configuration.  Arbi¬ 
trarily,  the  parameters  of  the  actual  configuration  correspond  to  the 
tabulated  values  of  the  parameters  with  the  exception  of  the  propeller 
plane  position.  Its  value  of  -0. 125  is  such  that  it  lies  midway  between 
two  tabulated  values. 

For  this  case,  three  evaluations  were  performed  with  x  =  -0.  25  , 

0.  0  and  0.  25  ,  the  corresponding  values  of  t  being  0.  3979  ft  , 

P 

0.4537  ft  ,  and  0.2969  ft  .  The  results  for  the  actual  configuration 
were  then  obtained  by  three -point  Lagrangian  interpolation  of  this  data. 
Figs.  6  and  7  are  the  completed  worksheets  for  x=  0.  25  . 

An  "exact"  evaluation  has  also  been  run  using  the  general  com¬ 
puter  programs.  The  two  results  are  presented  in  Fig.  8  along  with 
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SUB-PROCEDURE  1 

DETERMINATION  OP  PARAMETERS 
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GROSS  GCOMETKIC 
AND  rUCHT  I»A>t\METt.RS 
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SUB-PROCEDURE  3 

COMPUTATION  OF  SHROUD  SECTIONAL  RADIAL  FORCE  AND  MOMENT  COEFFICIENTS  AND  CENTER  OF  PRESSURE 


29-27x28 
d3«2lx32 
38  a  29*  30 

37  -  33  X  SS 

38  •  34  X  3& 

39  •  38  •  37  »  38 


COMPUTATION  OF  1 
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0.0815 

1 

•0.1115 

2 

-0.2200 
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Write  ihe  Krsi  three 
element*  of  coukx  81  as  the 
clement  of  couk'U  82.  83at¥J8f 
respectively. 


®  88^82x85 

89^83X87 

Cj^  •  89  •  88  ♦  88 


©  71  *83x70 

73-84x72 

c^*74-71.73 


Q  75- 74-^89 

Xcp/C  «  77  »  75  •  76 


82 

83 

84 
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0.297 

0.109 

85  88 

3.1418  »0.405 


87 

88 

89 

E^}— 1 

>1.5708 
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•0.872 
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71 
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73 

74 
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75 

78 
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Write  the  clement#  of  cotixN  5  a*  the  element*  of  coumv  79. 


W'rite  the  element*  of  rouw  33, 

correspcndini:  10  (^*-1.2.3...  11, 

#*  the  correspoodlnic  element* 
of  cotaw  83. 

©  Write  Ihc  elements  of  coumv33. 
CO  r  respond  init  to  i'-2,3,4...  12. 
as  the  corresponding  element# 
of  rouvs  88. 

Write  the  first  element  of 
coux*  33.  corresponding  to 
V  0  .  a*  the  element  of 
cix^xs  89. 


©  Write  the  elcm«  nt*  of 

cou-Mi  61.  corresponding 

toiriil,2.3  ...  II.  a*  the 
corrciponding  elements 
of  cou'Xx87. 

Write  the  first  element  of 
couvx  81.  corresponding  to 
v  >0  ,«*  the  element  of 
cou-xs  91, 

Write  the  second  element  of 
eou-xx  81  >  corresponding  to 
V‘‘\  .as  the  element  of 
couxs  92. 
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W’rile  the  second  element  of 
Cou-Mv  33,  corresponding  to 
-I  .  as  the  element  of 
couxx  90. 
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Write  the  element*  of  eouxsSI. 
corresponding  to  lf-2,3,4,..  12. 
at  the  corresponding  elements 
of  cou-xs  84. 
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98  «  90  X  92 
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FIGURE  6 

COMPLETED  WORKSHEET  I  FOR  INTERPOLATION/EXIRAPOI 


SUB-PROCEDURE  2 

COMPUTATION  OP  THREE-DIMENSIONAL  GUUERT  COEFOCIENTS  OF  EFFECTIVE  CAMBER 


Write  the  elemcAte  e(  cou-mkS  »•  the  (ir«t 
elemcoto(coU'UNi7.  9.  Ut  13,  IS,  17  and  19 
respectUcly, 


20««x7 

21  *9X9 

22  •  10X11 
22  •  12  X  13 
24  *14X13 
23*13X17 
23-13X19 

34  •  20  •  21  *  22  *  23  •  24  •  2&  *  23 


gaatiBrTBrMEBiE33WIMiBES2WB 
iB&^HDHDillllBQiiESSIllKHKBiSQi 


iPEUi 

Iksui 


Braann 

_ Mnwl 

KEEiX*!'!'Xj 

|^5B' 

■asil 


lEaHal 

IlSUJIj 

I 

IbemI 


_ — 

;EEll?tT3iBRTlBEl  D 

iB'B'HBPH'MBB.H.Wir'tW  Q 
- :BEa— BMP 


Vrito  the  lir«t  Acveft  clemeflts 
ol  cou-x-  39  as  the  first  element  of 
COUNHs  41.  43,  43.  47.  49.  31  and  33 
mpecti^ely. 
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38  *  48  X  49 
39-30x31 
60*32X53 


^^1 

iBCSal 

Bt!Wl 

BltHHlI 


PPPPBEESBPPBEESBUmBiniPPB^lPPE^CTl 

CEniBBBBS3IPBBB^UmnTTam4Mi*!MitiHi!wiB»aTaBn!!!!!s| 

llEI3lpBBlfi!;ii!2BpBkjLia*2!Slifflmii3i!23lBLli!iJI*’X'^^^*i^xa 


l»?.«.9...  19 


BBIQII 

BBSa 

nit™ 

deep' 

PEsail 

PEBM 

PEHOI 

PBESl 

PEQOi 

mpna 

liraB^gBl 


N/EXTRAPOLATION  EXAMPLE,  X  =  0. 25 
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FIGURE  7 


COMPLETED  WORKSHEET  H  FOR  INTERPOLATION/EXTRAP 
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QM  Write  the  ckmenu  of  coU'V'i  i$7  at  tbo 
correspondiDgelemenuof  cou‘X'>23S 
239  219*  220*  221  •  222  «  229  *  224  <  225 
240''  226  X  227 
242  >  237  X239 

Cp  tMItn  •  243  -  239  •  239  •  240  •  241  \  oUAnOAl 
cl  avM„  .  244 . 221  .  240 . 241  •  242 1  WADRAl 

Cptot-irw  -243 -239  *  239  *  239  •  2401 

-2«-2?3-?i?*24''»24»/  UNiJsn 


(7)  219^205x206 
220-207x208 

221  •  209  X  210 

222  -  211  X  212 

229'219x2I4 
224  .  215  x  216 
225*217X218 
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O  "EMct" 

□  /'cfuol, without  Propeller 


A  COMPARISON  OF  INTERPOLATED 
AND 

’’EXACT'’  SHROUD  SURFACE  PRESSURE  DISTRIBUTIONS 
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-0.30 

-0.20 

-0.10 

0 

0.10 

0.20 

0.30 

0.40 

0.50 


SHROUD  OUTER  SURFACE  PRESSURE  DISTRIBUTION 

Interpolated  "Exact 


x=-0. 125  x=-0. 125 


CO 

.241 
.107 
.063 
.017 
-0.039 
-0.068 
-0.071 
-0.008 
0.046 
0.132 


SHROUD  INNER  SURFACE  PRESSURE  DISTRIBUTION 

Interpolated  "Exact 


x=-0.125  x=-0. 125 


x=-0.25 

x=0.0 

x=0. 25 

-CO 

-oo 

-oo 

-1.051 
-0. 955 

-1.024 
-0. 955 

-0.837 

-0.931 

-1.075 

-0.760 

-0.876 

-1.039 

-0.704 

-0.746 

-0.983 

-0. 652 

-0.624 

-0.932 

-0.515 

-0.485 

-0.736 

-0.339 

-0.322 

-0.414 

-0.119 

-0.122 

-0.114 

0.154 

0.130 

0.220 

-  CO 
21 
39 
78 


-0.812 
-0.701 
-0.596 
-0.465 
-0.317 
-0. 122 
0.128 


-CO 

29 
55 


-0.913 

-0.781 

-0.673 

-0.618 

-0.485 

-0.326 

-0.120 

0.132 
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the  tabulation  of  the  three  sets  of  interpolation  data.  The  extent  of  the 
interpolation  involved  is  reflected  by  comparison  with  the  results  in¬ 
cluded  for  the  case  with  no  propeller  present.  Again,  we  find  that  the 
overall  agreement  is  very  good  and  cl  '  ’ly  seems  to  justify  such  inter¬ 
polation  or  extrapolation. 

2. 4  Example  for  Parametric  Studies 

The  third  and  final  example  illustrates  the  use  of  the  evaluation 
procedure  as  a  vehicle  for  parametric  studies.  We  have  chosen  to 
find  the  variation  of  the  shroud  thrust  coefficient  with  the  shroud 
incidence  holding  the  remaining  parameters  fixed. 

For  this  calculation,  it  is  not  necessary  to  consider  a  specific 
configuration  nor  SUB-PROCEDURES  3,  5  and  6  .  The  worksheets 
appropriately  completed  for  i  =  -0.20  ,  0.0  ,  and  0.20  rad  are 
shown  altogether  in  Fig.  9  .  Only  those  columns  having  more  than 
one  value  in  any  element  have  to  be  changed  in  performing  the  calcula¬ 
tion  for  each  value  of  i  .  The  uppermost  value  corresponds  to 
i  =  -0, 20  rad  ;  the  lower-left  value  corresponds  to  i  =  0. 0  rad  ; 
and  the  lower-right  value  corresponds  to  i  =  0.20  rad  . 

The  results  of  the  evaluation  are  plotted  in  Fig,  10  and  show  that 
the  variation  is  exactly  linear  over  the  range  examined.  This  result 
is  not  unexpected.  That  is,  the  shroud  thrust  coefficient  is  determined 


SUB-PROCEDURE  I 

DETERMINATION  OF  PARAMETERS 
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SUB-PROCEDURE  3 

COMPUTATION  OP  SHROUD  SECTIONAL  RADIAL  FORCE  AND  MOMENT  COEFnCIENTS  AND  CENTER  OF  PRESSURE 


«2  «3  f 


Write  the  Ural  three 
elcmceis  of  cou'X’i  61  aa  the 
element  of  eouvM  62,  63  and 
respective!)'. 


®  66*62x65 

68  ■  63  X  67 

e.«63-C6*68 


_  65  66 

•  I  -3.1416 


_  67  68  6» 

rT/Fl—  -I. $708 


(a)  7l*6:x70 

73  *  64  X  72 


©  75  •74*6? 

*c/C"77  *  75  *  76 


_  70  71 

ri/i  -0.3327 


— •  72  73  1  74 

t/8  0.3327  j  i 


75  76  77 

-0^ 


Write  the  elements  of  cou-vx  S  at  the  elements  of  coumx  73. 
Write  the  elements  of  cou-Mx  33.  Write  the  elements  of 


Write  the  elements  of  counx  33. 
corresponding  to  P'1,2,3...  11, 
as  the  corresponding  elements 
of  COU'MX  83. 

Write  the  elements  of  wawSS. 
corresponding  to  (''2.3.4  ...  12. 
as  the  corresponding  elements 

of  COU'MH  86. 

W’rite  the  first  element  of 
enux*  33.  corresponding  to 
1'  '0  .  as  the  clement  of 
couv-r  83. 

Write  the  second  element  of 
cotcMX  33,  corresponding  to 
P  '1  ,  as  the  element  of 
coals'*  30. 

Write  (he  elements  of  exMcxx  61. 
corresponding  to  ('•'2.3.4...  12. 

■!« 

of  eou-Mx  84. 


Write  the  elements  of 
cuuw  61.  corresponding 
(0C»1.2.3  ...  II.  as  the 
corresponding  elements 
of  couvx  87. 

Write  the  first  clement  of 
coui**  61.  corresponding  to 
(''*0  ,  as  the  element  of 
cou-M*  31. 

Write  the  second  clement  of 
cou'Mv  61.  corresponding  to 
p  •!  .  as  the  clement  of 
cou-Mv  32. 


80*78x73 
82*80  x81 
85  •  83  X  84 
88*  86x87 


®  34  *83X31 
^  96*  83X32 
38  •  30  X  32 

100  *  Algehrftic  sum  of  85 


FIGURE  9 

COMPLETED  WORKSHEET  I  FOR  PARAMET 
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SUB-PROCEDURE  2 

COMPUTATION  OP  THREE-DIMENSIONAL  OLAUERT  COEPHCIENTS  OP  CPPECTIVE  CAMBER 
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Write  the  element*  of  coummS  a*  the  (lr*t 
element  of  rou-Mts 7.  9,  ll.  13,  IS,  17  and  19 
respeelhely. 
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23  •  12  X  13 

24  •  14  X  1$ 

2$  •  16  X  17 
26  •  16  X  19 
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SUB-PROCEDURE  4 

COMPUTATION  OP  SHROUD  THRUST  COEFFICIENT 


S  a*  tlie  elements  of  coU'W'i  79. 

Write  the  element*  of 
couw  61,  correspondinic 
tot'll. 2, 3  ...  II.  a*  the 
correspondlAK  elements 
of  coct-xs  S7. 

Write  the  first  clement  of 
couv'i  61.  corresponding  to 
Vfi0,»*  the  element  of 
Cou-w*  91. 

W'rile  the  second  element  of 
cou-M*  61,  correspondlnKto 
i'  «1  ,  as  the  clement  of 
COUMf  92. 
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FIGURE  9 


FOR  PARAMETRIC  STUDY  EXAMPLE 
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by  the  product  of  the  three-dimensional  Glauert  coefficients  of  the 
effective  camber  and  the  two-dimensional  Glauert  coefficients  of  the 
propeller  induced  camber.  With  all  parameters  fixed,  except  i  ,  the 
three-dimensional  Glauert  coefficients  of  the  effective  camber  are 
independent  of  shroud  incidence.  Therefore,  the  shroud  thrust  co¬ 
efficient  will  also  be  proportional  to  the  shroud  incidence. 

The  crossover  point  on  Fig.  10  where  =  0*0  determines  the 
value  of  shroud  incidence  at  which  the  shroud  has  no  effect  on  the 
inflow  to  the  propeller.  For  angles  to  the  right  of  this  point  the  shroud 
will  increase  the  inflow  and  for  angles  to  the  left  of  this  point  it  will 
decrease  the  inflow  to  the  propeller. 
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SUB-PROCEDURE  1 
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DETERMINATION  OF  PARAMETERS 
(Use  WORKSHEET  I) 

1.  For  the  configuration  to  be  evaluated,  enter  N  ,  the  NACA  airfoil 
section  designation,  and  the  NACA  mean  line  designation  in  the 
appropriate  locations  below  the  figure  on  WORKSHEET  I  .  Also 
fill  in  all  the  dimensions  c,  R,  . . .  ,  indicated  on  the  figure  it¬ 
self  with  these  e>r.ceptions: 

i.  and  are  applicable  only  when  an  NACA  4-digit 
mean  line  is  being  used. 

ii.  and  the  inset  in  the  figure  on  WORKSHEET  I  are 
applicable  only  when  an  NACA  6-digit  mean  line  is 
being  used. 

The  following  sign  conventions  must  be  observed: 

i.  Xp  ,  x^  and  x^  are  positive  or  negative  depending 
on  their  location  relative  to  the  chosen  coordinate 
system,  i.  e. ,  positive  to  the  right  of  the  midchord 
and  negative  to  the  left. 

ii.  y_  is  measured  relative  to  the  chord  line,  and  is 
positive  when  above  and  negative  when  below  this  reference. 

iii.  All  other  dimensions  are  positive. 
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iv.  i  is  positive  in  a  counter-clockwise  sense. 

Record  the  physical  flight  conditions  T  ,  U  ,  and  p  . 

2.  Compute  the  remaining  parameters  from  the  equations  which  are 
given  and  the  data  from  step  1.  Write  the  results  in  the  columns 
headed  actual.  The  following  exceptions  should  be  noted: 

i.  p  and  m  are  applicable  only  when  an  NACA  4 -digit 
mean  line  is  being  used. 

ii.  a  and  c ^  are  applicable  only  when  an  NACA  6-digit 

i 

line  is  being  used,  c^  is  presumably  known  and  is 

i 

simply  written  down;  it  is  positive  if  the  mean  line  is 
concave  downward  (^)  and  negative,  if  concave  upward 

. 

3.  Use  the  TABLE  OF  PARAMETERS  and  choose  where  necessary 
the  values  of  the  tabulated  parameters  which  are  closest  to 
(larger  or  smaller)  or  equal  to  the  actual  values;  write  these 
values  in  the  indicated  locations  under  the  column  headed 
EVALUATED.  For  the  parameters  which  can  have  arbitrary  values, 
use  the  actual  values  in  the  evaluation  except  for  t^/c.  The 
value  to  use  for  t^/c  is  described  in  the  next  step. 

4.  For  finite  thickness,  compute  the  consistent  value  of  t^/c  to  be 
used  in  the  evaluation  as  follows: 
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i.  Use  the  value  of  e  ,  pt  and  X  from  the  column  headed 
EVALUATED  ,  and  find  the  corresponding  shroud  thickness 
ratio  at  the  propeller  plane  t  /c  from  TABLE  2. 1  . 

ii.  Multiply  the  value  of  x  in  the  column  headed  evaluated 
by  the  value  of  c  and  obtain  a  "new"  value  of  x 

p 

iii.  Using  the  actual  configuration  determine  the  "new"  value 
of  tp  corresponding  to  the  "new"  value  of  . 

iv.  Take  t  from  the  figure  on  WORKSHEET  I  and  with  the 

"new"  value  of  t  form  the  ratio  t  /t  and  multiply  it 
p  m  p  ^ 

by  the  value  of  t  /c  from  TABLE  2. 1  .  Enter  the  result 
in  the  location  t^/c  in  the  column  headed  evaluated  . 

For  zero  thickness,  use  the  value  0  in  the  evaluation. 

5,  In  the  remaining  sub-procedures,  the  values  of  all  parameters  re¬ 
quired  will  be  taken  from  the  column  headed  evaluated. 
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SUB-PROCEDURE  2 

COMPUTATION  OF  THREE-DIMENSIONAL  GLAUERT  COEFFICIENTS 

OF 

EFFECTIVE  CAMBER 
(Use  WORKSHEET  I) 

1.  Write  the  value  of  in  column  2  . 

2.  Write  the  value  of  C,p  in  column  32  . 

3.  Write  the  appropriate  geometric  camber  parameters  in  the 
indicated  columns  as  follows: 

i.  For  an  NACA  4-digit  mean  line  write  the  value  of  m  in 
COLUMN  28  and  the  value  2i  in  column  30  . 

ii.  For  an  NACA  5-digit  mean  line  write  the  value  1  in 
COLUMN  28  if  the  mean  line  is  concave  downward,  or  the 
value  -1  if  the  mean  line  is  concave  upward;  also  write 
the  value  2i  in  column  30  . 

iii.  For  an  NACA  6-digit  mean  line,  write  the  value  C£  in 

i 

COLUMN  28  and  the  value  2i  in  column  30  . 

4.  Use  the  value  of  A  nnd  the  thickness  parameters  to  pick  the 
coefficients  of  the  shroud  source  strength  distribution  a^  from 
TABLES  3. 1  -  3. 2  and  proceed  as  follows: 

i.  For  an  NACA  4-  or  5- digit  airfoil,  write  the  appropriate 
column  of  numbers  from.  TABLE  3. 1  in  column  1  and  the 


value  1  in  column  3  . 
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ii.  For  an  NACA  6- digit  airfoil,  choose  the  column  headings 

in  TABLE  3. 2  corresponding  to  the  6 -digit  family  with 

the  same  or  closest  smaller  value  of  t  /c  and  write 
-  m 

this  column  of  numbers  in  column  1;  if  the  table  does  not 

contain  a  smaller  value  of  t  /c  ,  use  the  column  with 

m  ’ 

the  smallest  value  of  t  /c  available.  Also  write  the 

m 

value  of  t^/c  from  the  column  heading  used  in 
COLUMN  3  . 

iii.  For  zero  thickness,  write  the  value  0  in  all  the  elements 
of  COLUMN  34  . 

5.  Use  the  val,ue  of  X  to  pick  the  elements  for  the  evaluation  of  the 
camber  induced  by  shroud  thickness  ^  from  TABLES  4. 1  - 
4.  3  .  Write  each  column  of  numbers  from  the  appropriate  table 
in  COLUMNS  6,  8,  10,  12,  14,  16  and  18  respectively.  For  zero 
thickness,  omit  this  step. 

6.  Use  the  geometric  camber  parameters  to  pick  the  two- 
dimensional  Glauert  coefficients  of  the  geometric  camber  from 
TABLES  5. 1  -  5.  3  .  Write  the  column  of  numbers  from  the  ap¬ 
propriate  table  in  column  27.  For  NACA  4-digit  mean  lines,  with 

p  =  0. 40,  0.  30,  0.  25,  0. 20,  or  0. 10  use  the  data  for  the  corre¬ 
sponding  negative  value  of  p  with  the  sign  of  the  numbers  for 
p  =  0,  2,  4,  6,  8,  10  and  12  reversed. 


7.  Use  the  values  of  the  gross  geometric  and  flight  parameters  to 
pick  the  two-dimensional  Glauert  coefficients  of  the  propeller  in¬ 
duced  camber  e„,  from  TABLES  6. 1  -  6. 18  .  Write  the 

F'  u 

column  of  numbers  from  the  appropriate  table  in  column  31  .  For 
X  =  0. 25  ,  use  the  tables  for  x  =  -0.  25  with  the  sign  of  the  num¬ 
bers  corresponding  to  u  =  1,  3,  5,  7,  9  and  11  reversed. 

8.  Use  the  value  of  X  to  pick  the  elements  for  the  curvature  cor¬ 
rection  of  thr  two-dimensional  Glauert  coefficients  of  the  effective 


shroud  camber  Oj^  ^  from  TABLES  7. 1  -  7.  3  .  Write  each  col¬ 
umn  of  nup'bers  from  the  appropriate  table  in  columns  40,  42,  44, 
46,  48,  50  and  52  respectively. 

9.  Compute  the  three-dimensional  Glauert  coefficients  of  the  effective 
3D 

camber  c^  by  following  the  instructions  given  in  SUB¬ 
PROCEDURE  2  on  WORKSHEET  I;  for  zero  thickness,  omit  the 
instructions  in  steps  1,  2  and  3  . 
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SUB-PROCEDURE  3 
COMPUTATION  OF 

SHROUD  SECTIONAL  RADIAL  FORCE  AND  MOMENT  COEFFICIENTS 

AND 

CENTER  OF  PRESSURE 
(Use  WORKSHEET  I) 


1.  Compute  the  shroud  sectional  radial  force  coefficient  c  ^  ,  the 

pitching  moment  coefficient  c  ^  ,  and  the  ncn-dimensional  center 

of  pressure  x  /c  by  following  the  instructions  given  in  SUB- 
cp 

PROCEDURE  3  on  WORKSHEET  I  . 


SUB-PROCEDURE  4 


COMPUTATION  OF  SHROUD  THRUST  COEFFICIENT 
(Use  WORKSHEET  I) 


1.  Use  the  values  of  the  gross  geometric  and  flight  parameters  to  pick 

the  terms  of  the  propeller-shroud  thickness  contribution  to  the 
shroud  thrust  coefficient  from  TABLES  8. 1  -  8.  27.  Write 

V 

the  column  of  numbers  from  the  appropriate  table  in  column  78  . 

For  zero  thickness,  omit  this  step  and  write  the  value  0  in  all  the 
elements  of  column  82  . 

2.  Write  the  value  of  C,p  in  column  81,  the  value  of  [i  in  column  111, 
and  the  value  of  in  column  112  . 

3.  Compute  the  shroud  thrust  coefficient  C^  by  following  the  in¬ 
structions  given  in  SUB-PROCEDURE  4  on  WORKSHEET  II  ;  for 
zero  thickness,  omit  the  instructions  in  step  1  . 
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SUB-PROCEDURE  5 

COMPUTATION  OF  NET  SHROUD  PRESSURE  COEFFICIENT 

(Use  WORKSHEET  H) 

Compute  the  net  shroud  pressure  coefficient  ACp(x/c)  by  follow¬ 
ing  the  instructions  given  in  SUB-PROCEDURE  5  on  WORKSHEET  II  . 
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SUB-PROCEDURE  6 
COMPUTATION  OF 

OUTER  AND  INNER  SHROUD  SURFACE  PRESSURE  COEFFICIENTS 

(Use  WORKSHEET  II) 

1.  Write  the  value  of  C^  in  column  227  . 

2.  Write  the  value  of  t^/c  in  column  228  . 

3.  Use  the  value  of  x  to  pick  the  elements  for  the  shroud  curvature 
contribution  to  the  shroud  surface  pressure  arising  from  the 
effective  shroud  camber  Sj^  ^  from  TABLES  9. 1  -  9.  3  .  Write 
each  column  of  numbers  from  the  appropriate  table  in  columns  160, 
162,  164,  166,  168,  170  and  172  respectively. 

4.  Use  the  value  of  x  to  pick  the  elements  for  the  shroud  curvature 
contribution  to  the  shroud  surface  pressure  arising  from  the 
shroud  thickness  Tj^  ^  from  TABLES  10. 1  -  10.  3  .  Write  each 
column  of  numbers  from  the  appropriate  table  in  columns  174,  176, 
178,  180,  182,  184  and  186  respectively.  For  zero  thickness, 
omit  this  step. 

5.  Use  the  values  of  the  gross  geometric  and  flight  parameters  to 
pick  the  direct  propeller  contribution  to  the  shroud  surface  pres¬ 
sure  coefficients  c^,(x/c)  from  TABLES  11. 1  -  11. 18  .  Write 
the  column  of  numbers  from  the  appropriate  table  in  column  226. 


^2 


For  X  =  0.  25  ,  use  the  tables  for  x  =  -0.  25  invert  each 

column  and  reverse  the  sign  of  each  number. 

Use  the  values  of  the  thickness  parameters  to  pick  the 

linearized,  two-dimensional  contribution  of  shroud  thickness 

2D 

to  the  shroud  surface  pressure  coefficient  c^  (x/c)  from 
TABLES  12. 1-12.2  as  follows: 

i.  For  an  NACA  4-  or  5-digit  airfoil,  choose  the  column 
heading  in  TABLE  12. 1  having  the  same  or  closest  smaller 
value  of  t^/c  and  write  this  column  of  numbers  in 
COLUMN  231  .  Also  write  the  value  of  tj^/c  from  the 
column  heading  used  in  column  229  . 

ii.  For  an  NACA  6-digit  airfoil,  choose  the  column  headings 

in  TABLE  12.  2  corresponding  to  the  6-digit  family  with 

the  same  or  closest  smaller  value  of  t  /c  and  write  this 
-  m 

column  of  numbers  in  column  231  ;  if  the  tables  do  not  con¬ 
tain  a  smaller  value  of  tj^/c  ,  use  the  column  with  the 
smallest  value  of  t^/c  available.  Also  write  the  value 
of  tj^/c  from  the  column  heading  used  in  column  229  . 

iii.  For  zero  thickness,  write  the  value  0  in  all  the  elements 
of  columns  195,  196,  197,  198,  199,  200  and  201. 

Compute  the  net  shroud  surface  pressure  coefficient  by  using 
SUB-PROCEDURE  5  . 
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8.  Compute  the  shroud  outer  and  inner  surface  pressure  coefficients 

c  (outer)  and  c  (inner)  respectively  by  following  the  instruc- 
P  P 

tions  given  in  SUB-PROCEDURE  6  on  WORKSHEET  H  ;  for  zero 
thickness,  omit  the  instructions  in  steps  2,  4  and  8.  In  this  sub¬ 
procedure  the  quadratic  form  of  the  two-dimensional  thickness 
contribution  to  the  shroud  surface  pressure  is  used.  If  it  is  de¬ 
sired  to  use  the  linearized  form,  then  follow  the  last  instruction 
of  STEP  9  . 


TABLE  OF  PARAMETERS 


NACA  5-Digit  Mean  Line 


Accuracy:  ±  0.  001 


5-DIGIT 

MEAN 

LINES 


TWO-DIMENSIONAL  GLAUERT  COEFFICIENTS  OF  PROPELLER  INDUCED  CAMBER 
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